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Abstract 

White dwarf stars composed of carbon, oxygen or heavier elements are 
expected to crystallize as they cool down below certain temperatures. Yet, 
simple arguments suggest that the helium white dwarf cores may not solidify, 
mostly because of zero-point oscillations of the helium ions that would dissolve 
the crystalline structure. We argue that the interior of the helium dwarfs 
may instead form a macroscopic quantum state in which the charged helium- 
4 nuclei are in a Bose-Einstein condensate, while the relativistic electrons 
form a neutralizing degenerate Fermi liquid. We discuss the electric charge 
screening, and the spectrum of this substance, showing that the bosonic long- 
wavelength fluctuations exhibit a mass gap. Hence, there is a suppression at 
low temperatures of the boson contribution to the specific heat - the latter 
being dominated by the specific heat of the electrons near the Fermi surface. 
This state of matter may have observational signatures. 



1 Introduction 



Consider a system of a large number of neutral atoms, each having Z electrons, and 
for simplicity we focus on Z < 10. Suppose we increase the number- density of the 
atoms - denoted by Jo/Z, where J is the electron number-density - so that the 
average separation between the atomic nuclei 
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becomes much smaller than the Bohr radius, but is still much greater than the 
size of the charged nuclei. This would certainly be the case when, for instance, 
J ~ (0.5 — 5 MeV) 3 , for which the electron Fermi momentum, pp = (Sn 2 ^) 1 ^ 3 , 
ranges from well-relativistic to ultra-relativistic. 

Under these circumstances the electron Fermi energy Ep will exceed the atomic 
binding energy, and the neutral atoms will dissolve into the electrons and nuclei. 
What is the ground state of such a system, when it's held together by gravity? 

The answer would depend on its temperature T. At certain high temperatures a 
plasma of negatively charged electrons and positively charged nuclei (ions) would be 
the equilibrium state. However, for lower temperatures the system would undergo 
significant charges. 

Let us consider the cooling process by ignoring the quantum mechanical effects 
for the nuclei (ions), while treating the electrons as a degenerate Fermi liquid^. 
Consider the nuclei that are bosons. The strength of interactions of a classical 
plasma of the charged nuclei (ions) is customarily characterized by the ratio of the 
average Coulomb energy of a pair of ions to the thermal energy (see, e.g., [T]) 

r = Ecouiomb _ {Zef /(A^d) ^ 

^Ephermal/^ kp;T 

here e denotes the electric charge, kg is the Boltzmann constant, and the second 
equality in assumes the validity of the classical approximation. 

It has been known that when T > 180 the ion plasma is coupled strongly- 
enough for the system to crystallize, i.e., when the temperatures become low enough 
and T reaches the above-mentioned value, the system would undergo crystallization 
transition (for earlier works, see Refs. [2j[3], for latter studies, [UES E] and references 
therein). This has direct relevance to white dwarf stars (WDs): It is expected that 
in most of the white dwarfs, consisting of carbon, oxygen, or heavier elements, the 
crystallization transition takes place in the process of cooling (for a review, see, e.g., 
Ref. 0). 

The above discussions were classical. Although it was observed in Ref. [8] 
that the quantum effects of the ions may be important for cooling of the WDs, 
quantitative studies using a model for a quantum solid and Lindemann's empirical 
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rule lead to a conclusion [6] that the melting curves for most of the WDs are not 
changed significantly by the quantum effects. 

In the present work we're interested in the system of electrons and helium-4 
nuclei. The latter being lighter than the other nuclei, and having lower charge, 
should be expected to reveal stronger sensitivity to the quantum effects. 

Such a system should be present in cores of helium white dwarfs. These dwarf 
stars stayed relatively cool due to their binary companions, which kept accreting a 
part of the dwarf's energy during their paired evolution. As a result, the temperature 
in such dwarfs was never enough to begin the helium burning process. 

The electron number-densities in the interval Jo — (0.5 — 5 MeV) 3 would cor- 
respond to the mass density of the helium-electron system 5 • (10 7 — 10 10 ) g/cm 3 . 
The densities closer to the lower end of this interval (which we refer to as the low 
densities hereafter) may well exist in cores of the helium white dwarfs. We will also 
discuss the densities closer to the upper end, as they may be present in some other 
astrophysical circumstances. 

The system should be long lived as the helium-4 nuclei are stable w.r.t. fission. 
Furthermore, some nuclear reactions that could contaminate the helium-4 cores by 
their products are suppressed. One of this is the neutronization process due to 
the inverse beta-decay. In our case the electrons with Jo ~ (0.5 — 5 MeV) 3 are 
not energetic enough to reach the neutronization threshold of the helium-4 nucleus, 
which is about 20 MeV. Moreover, we would expect that the rate of helium fusion via 
the triple alpha-particle reaction is suppressed as the so-called pycnonuclear reaction 
(i.e., nuclear fusion due to zero-point oscillations in a high density environment) 
rates are exponentially small [TJ. Hence, the cores of these WD's are expected to be 
dominated by helium-4 for very long time after their formation. 

In what follows, we will focus on the case when Z = 2. Then, for Jo ~ (0.5 — 
5 MeV) 3 , the classical considerations of Eq. (T5]) would give the crystallization 
temperature of the helium ions T cryst ~ (10 6 — 10 7 ) K. However, before these 
temperatures are reached from above, the de Broglie wavelengths of the helium 
nuclei would start to overlap, suggesting that the quantum effects of the ions may 
be important. An estimate for the critical temperature below which this will happen 
can be given by 

4tt 2 

For the number densities at hand we obtain T c ~ (10 7 — 10 9 ) K. These are greater 
than the crystallization temperatures estimated above. For this, it's unjustified to 
approximate the thermal energy of the ions by its classical expression 3k bT 72, as it 
was done in the denominator of ([2]n. 



2 Note that for the system of carbon nuclei we would obtain the crystallization temperature 
T cryst ^ 6 • (10 6 - 10 7 ) K, and the corresponding critical temperatures T c ~ 2 ■ (10 6 - 10 s ) K. 
Thus, for the low densities the crystallization temperature exceeds the critical one, and the carbon 
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Perhaps the most important quantum effect, that invalidates the predictions of 
d2J), is related to the zero-point fluctuations: Even at very low temperatures, close to 
the absolute zero, there will be a significant amount of energy stored in the quantum 
mechanical zero-point oscillations of the nuclei in a would-be crystal. This energy 
can be approximated as follows (see, e.g., pQ): 

* = ^ s §e£ = (3^)% 3 ao'-io').v, (4) 

where tuh — 3.7 GeV is the helium nucleus mass, and u l ° ns denotes the ion plasma 
frequency. This energy is commensurate with the classical thermal energy at tem- 
peratures (2/3) (3 • 10 7 — 10 9 ) K - the values that are greater than the crystallization 
temperatures obtained from (j2j). 

The above discussions suggest that a measure of coupling in such a cold plasma 
is the ratio of the Coulomb energy to the zero-point energy. This ratio, properly 
normalized to be consistent with the conventions of [lj (which are different from 
ours), takes a small value that is not enough for the plasma to crystallize. It increases 
with the nuclear charge Z, and the atomic number A. Thus, for heavier nuclei, it 
will be high-enough not to obstruct the crystallization. In the charged helium-4 
case with the densities at hand, however, the zero-point oscillations would melt the 
crystal latticdjl 

The above considerations lead us to conclude that the dense system of the helium- 
4 nuclei and electrons may not solidify, mostly because of the zero-point oscillations 
of the helium nuclei. 

2 Charged condensation 

What could then be an adequate description of the quantum state discussed above? 
It was argued in Ref. JO] that a state in which the charged helium-4 nuclei form 
a Bose-Einstein (BE) condensate, while the degenerate electrons neutralize the net 
charge of the condensate, represents a local equilibrium state of the system. 

solidification would take place. For the higher densities the opposite statement would be true. 
One should note, however, that we're comparing these for a fixed electron number-density in the 
helium and carbon cases, which implies three times greater mass density for the carbon case, so 
the comparisons should be readjusted by this factor. Since the carbon cores is not a subject of the 
present paper, we will not deal with this procedure here. 

3 The case for the crystallization is worsened further by the fact that for the relativistic electrons 
with the Fermi energy Ep ~ (37r 2 Jo) 1 ^ 3 , the Coulomb interaction gets screened. The screening in 
such a plasma becomes effective at scales of the order the Debye wavelength, \r> ~ (Ep / 'e 2 ''Jo) 1 / 2 ~ 
(e 3 J r 1/3 , and can be appreciable at the scale of the average inter-particle separation d. Hence, 
not only the denominator in ((2|) should be modified, but also the numerator should be reduced 
(multiplied by a factor) to account for the screening. The question of how significantly in the 
(ultra)relativistic case this screening may affect the crystallization of the carbon nuclei is left open. 
Note that this screening can be important in a nonrelativistic setup as well, see, e.g., Ref. g]. 
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The helium nuclei could have condensed to such a state when their de Broglie 
wavelength started to overlap, which according to our estimates above happened 
before the temperature dropped to the would-be crystallization value. Such a state 
of the helium nuclei cannot be resolved into individual particles. It represents a 
macroscopic mode with a large occupation number, somewhat similar to a quantum 
liquid of charged particles, where the charge is screened. This state, referred to as 
charged condensate, has a number of distinctive signatures: a photon propagating 
through this medium slows down, since it acquires a Lorentz-violating mass term, 
and its longitudinal polarization exhibits unusual dispersion relation. Moreover, 
electric charges in this medium, were shown in [9], are screened at the distance 
scales 1/M, where M for the helium-electron system denotes 

M = (2e 2 m H J ) 1/ \ (5) 

Importantly, the screening takes place at scales larger than 1/M; on the other hand, 
1/M < d (see discussions of this in Section 3)! 

In Ref. [9], however, the dynamics of the electrons was neglected, since the 
fermions considered there were heavier. One of the goals of the present work is 
to include the dynamics of electrons. We will show that the screening length can 

1 /3 

be approximated by as long as m# ^> J (as it is the case here). However, 
we will also see that, inclusion of the electron dynamics modifies quantitatively the 
spectrum of small perturbations above the ground state^. 

Before turning to these issues, let us deal with the gravitational part of the 
problem. From now on we consider the system at temperatures much lower than ~ 
10 7 K, and ignore thermal effects. Thus, we discuss the electrons in the gravitational 
and electrostatic fields in the interior of an astrophysical object, e.g., in the core of 
a helium dwarf. The Fermi momentum, pp, is related to the chemical potential /i 
as follows: 

\]pf + K + V g + eA = fi, (6) 

where m e is the electron mass, V g is the (negative) gravitational potential energy 
per electron, and A denotes the gauge potential. An average A Q would be nonzero 
for a nonzero net surface charge. However, we will assume that the net charge, 
even if present, is small enough, so that \eAo\ <C \V g \. Then, as is well known, the 
overall stability of a dwarf star is due to the balance between attractive gravity and 
repulsive Fermi degeneracy pressure of the electrons. Thus, we briefly discuss the 
interplay between V g and the Fermi energy. 

4 It was argued in Ref. [10], that there could exist long-lived spherically symmetric objects of 
a (sub) atomic size with the charged condensate of helium-4 nuclei and electrons in their interior 
that could be held together by electrodynamic forces, as long as there is a significant amount of 
uncompensated charge at the surface of such objects. This interesting possibility and it applications 
will be studied further elsewhere. Here, the force that keeps the system together is gravity. 
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For this one could use the Thomas- Fermi (TF) approximation: Ignore A in ([6]), 
and deduce from it pp = {{ji — Vg) 2 —m 2 e ) 1 / 2 . Relate the latter to the number density 

and use this expression in the Poisson equation for V g to derive the Chandrasekhar 
equation^. 

Our main interest is in local properties of particles at scales that are much 
smaller than the size of the star. These properties are determined by electromag- 
netic interactions, once the average effective chemical potential for the particles in 
a given homogeneous region of the star is set by the gravitational dynamics. For 
instance, consider a small region in the core of the star. Local distortions in the 
number-density (i.e., distortions in the local Fermi momentum pp(x)) will cause the 
respective change in the local field potential Aq(x), while the change in the local 
gravitational field is negligible since the gravity is so much weaker. 

These arguments can be formalized by the following relation: 



\/p 2 f{ x ) 



+ ml + eAo(x) = (fi - V g (x)) = fif , (8) 



where the brackets, (...), denote the averaging of the large-scale gravity effects, 
over the region of a roughly uniform number density in the star. Hence, by fif we 
denote the effective chemical potential that governs the fermion dynamics in a local 
approximately homogeneous region, e.g., in the core of the star. The size of the 
homogeneous region should be greater than the photon Compton wavelength, for 
our discussions below to be valid. 

Having the gravity part dealt with, let us turn to the microscopic non-gravitational 
physics. For a large number of quanta that form a macroscopic state, such as the 
charged condensate, the description in terms of the effective field theory of the order 
parameter, and its long wavelength fluctuations, becomes adequate. The effective 
field theory will be constructed based on the fundamental principles of symmetries 
and properties of the interactions involved. This is the formalism that we will adopt 
for the rest of the work. 

The effective Lagrangian that captures the microscopic properties of the system 
contains the photon field A^, a charged scalar field 0, and fermions with the 

effective chemical potential [if discussed above 

£ = ~\ F lu + I A^l 2 - + t(i7"£>„ - m e )^ + ^tt+tt . (9) 

In a conventional setup, upon quantization, the fluctuations of <fi would describe 
the helium-4 nucleus of charge +2e, and its antiparticle. Such a description would 



5 Since this equation and its consequences are well studied, we'll not discuss them here. For its 
derivation in the TF approach, and further details see, Ref. and references therein. 
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be valid at densities that are above the atomic, but well below the nuclear scales. 
Here, instead, the vacuum expectation value of the field <fi will serve as an order 
parameter for the condensation of the helium-4 nuclei, thus describing a state with 
a large occupation number. Because of this large number, the description of such 
a quantum state in terms of the classical field is valid. Fluctuations of the order 
parameter will be discussed below. 

The fermions in (JH]) describe the (quasi) electrons with charge — e, and (quasi) 
positrons. The covariant derivative for them is defined as _D M = <9 M + ieA^, while for 
the scalar it reads 

Dp = — 2ieA tM = d„- i{2eA tl + fij^) . (10) 

The latter is equivalent to the introduction of the chemical potential fi s for the 
helium-4 nuclei. 

The Lagrangian has a global U s (l) symmetry, yielding the conservation of 
the number of scalars, as well as another global Uf(l), responsible for the fermion 
number conservation. One linear combination of these two symmetries is gauged, 
and the corresponding conserved current is coupled to the photon field. We could 
add the quartic scalar self- interaction term to (Q, but this won't change our results 

1/3 

significantly, as long as the quartic coupling is not strong, and m# 3> J . 

Since the electrons are relativistic, the quantum loops would renormalize the 
effective Lagrangian by generating additional terms proportional to the fermion 
chemical potential and temperature, which would be still consistent with the sym- 
metries of the theory. These terms would come suppressed by the fine structure 
constant, e 2 /4ii, and we do not expect them to modify our quantitative results sig- 

1/3 

nificantly in the regime when m# ^ J . Hence, for simplicity of presentation we 
do not include those terms here. 

For further convenience, we introduce the following notation for the scalar field, 
= ^c"e 4a , and work in the unitary gauge where the phase of the scalar is set to 
zero, a = 0. In this gauge, the Lagrangian density reads: 

£ = ~F% + \{d,af + { ^Ay - l -mW + VirfDp - m e )* + . (11) 

The kinetic term for the scalar in (Q gives rise to the third term in the Lagrangian 
(|lip . If the field A Q acquires an expectation value, this term plays the role of a 
tachyonic mass for the scalars (see, e.g., [12]). In particular, when (2eA ) = rnn, 
the scalar field condenses. That is to say, the equations of motion that follow from 
( TTTj) have the following static solution: 



{2eA ) + n 8 = m H , ( a ) = J . (12) 



2m 



H 



What fixes the value of Aq is the net surface charge. In particular, when the surface 
charge is absent, (2eA ) = 0, and has to be /i s = run, for the condensation to occur. 
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The system is neutral in the bulk, as the charge in the condensate is exactly canceled 
by the electron charge |9j. 

Dynamically, the condensation proceeds in a conventional manner: At tempera- 
tures higher than T c most of the states are in thermal modes, and n s {T) is less than 
run, but increases as the temperature dropa^j. As T — > T c a significant fraction of 
the modes ends up in the zero momentum ground state, while /i s (T c ) asymptotes to 
tuh (for a discussion of relativistic BE condensation, see, e.g., [T5]). 

The last comment in this section is that we're looking at the effects of the con- 
densate that are l/m H suppressed. In particular, the order parameter of the con- 
densation that we're discussing (given in eq. ({12]) ) is the l/m H suppressed quantity. 
Clearly, in order for the condensate to make sense at some fixed nonzero temperature 
T, the de Broglie wavelength of the bosons, ~ l/Vm^T, should be larger than the 
inter-particle separation - the condition which would break down for very large m# 
(in agreement with the decoupling), but is fulfilled for the values of the parameters 
in our case. Furthermore, the inverse of the order parameter (a) should be smaller 
that the physical size of the system - the condition that would also be violated for 
some very large ra# (consistent with the decoupling), but is preserved in our case. 



3 Charge screening and small perturbations 

As a next step we'd like to discuss screening of a probe electric charge in the conden- 
sate. To determine the screening length we consider a small, spherically symmetric 
object with a nonzero charge placed in the condensate. Outside of the charge, in 
the condensate, the equations of motion for a static A and a, as derived from (J9]) 
are: 

- V 2 i + 4eVi = e J , - VV = (4e 2 i 2 - m 2 H )a . (13) 
Following [9] we parametrize the perturbations as follows: 



Mr) = |f + SA (r) , a{r) = + 5a{r) . (14) 

To include the effects of the fermion fluctuations we substitute Jo — (Pf/37t 2 ) on the 
r.h.s. of the first equation in (fT3"|) . and express pp in terms of the gauge potential 
A Q , and the effective chemical potential [if via ([8]). The result is that 5Jo(r) gets 
related to 5Aq. Thus, the coefficient in front of SAq in the final equation for the 
perturbations gets modified as compared to [9] 



- V 2 5A + e 



2J | (37t 2 Jq) 2 / 3 
run vr 2 



5A = -2M 2 Sa, (15) 



6 Here we ignore the temperature dependence of tuh which would be present due to the quantum 
loop effects. 
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-V 2 5a = 2M 2 5A . (16) 

The above approximations is valid as long as we focus on solutions that satisfy 
5a <C a/ Jo/2rriH and 5A <C tuh /2e. 

In the regime when mjj» Jq^j we find that M 2 ^> J^ 3 , and then we can neglect 
the second term on the l.h.s. in (115p . For large r we require that <L4q, <5cr — > 0. The 
solutions are similar to those of [9], in which the boundary conditions select the 
decaying functions: 

e -Mr 

5A (r) = — — [ci sin(Mr) + c 2 cos(Mr)] , (17) 

Q -Mr 

Sa(r) = — — [-ci cos(Mr) + c 2 sin(Mr)] . (18) 

The constants c\ and c 2 are to be determined by matching these solutions to those 
in the interior of the small, charged object. Thus, for a probe particle, the screening 
occurs at scales greater than 1/M. For the helium-electron system at hand, 1/M 
is larger than the size of the helium nuclei, but is shorter that the average inter- 
particle separation d. The latter statement needs some qualifications. The effective 
field theory description adopted above is not expected to hold at scales below d. 
There is an explanation of this in terms of a cancellation between the potentials 
due to the two long-wavelength modes - Coulomb and "phonon" quasiparticles - 
both of which are much lighter than 1/d, and are well-within the validity of the 
effective field theory. This and related issues will be discussed in detail elsewhere. 
Here we only consider phenomena at scales greater than d, where screening (fTTj) . 
(fl8"j) is significant, and the effective field theory descriptions is valid. 

The magnetic interactions for the fermions filling the Fermi sphere, are not 
screened at scale M~ x . Instead, as it will be clear from the Lagrangian (1191) derived 
below, the magnetic interactions are screened at scale of the Compton wavelength 
of the massive photon. This scale will end up being much greater that the average 
inter-particle separation d. 

Let us now turn to the spectrum of small perturbations. The electrons fill the 
Fermi surface with the Fermi energy Ep. As to the bosons, there are three com- 
ponents of the massive photon, and the heavy scalar mode. We "integrate out" 
the fermions by taking into account their effect via the TF approximation, like we 
did it above for the static case. The resulting Lagrangian density for the bosonic 
perturbations takes the form: 

Ax, = ~\fl v + \ml{5A Q f - ^m 2 (5AA 2 + \{d,5af + 2M 2 5A 5a, (19) 
where f^ u denotes the field strength for dA^, j = 1, 2, 3 is the spatial index and 

e 2 (3vr 2 J ) 2 / 3 ^ 



m 2 = m 2 + ^ ^-jp , m 7 = 2eW^-. (20) 
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The calculation of the spectrum is a bit lengthy but straightforward. The result is 
as follows: The two (out of three) helicities of the massive photon are transverse 
and propagate according to the massive dispersion relation 

uj 2 = k 2 + m 2 , (21) 

where k denotes the spatial three momentum. Furthermore, there are two modes 
with the dispersion relations 

9 ,of m l + m l\ 2M 4 ml 

ujI = k 2 — - — =-1 H =- H 1 

± V 2m 2 , y ml 2 

± W4k2 I l + 

V m o 

The solution with the plus subscript corresponds to the scalar mode and its mass 
squared in this frame is a>+(k = 0) = (4M 4 /7tt,q), while the solution with the minus 
subscript corresponds to the longitudinal component of the massive vector field, 
with cu 2 (k = 0) = m 2 . Here it is useful to note that the expression for M (jSJ) can 
be rewritten as M = (m^m 7 ) 1 / 2 . 

Furthermore, in the limit when the fermions are non-dynamical (i.e., are frozen 
"by hand" or some other dynamics), then m — > m 7 , and the solutions reduce to the 
ones obtained in [5J, however, in the present setup the difference between m and m 7 
is greater than m 7 . Therefore, the fermion dynamics contributes by an additional 
screening of the electrostatic interactions. 

The solutions (T2~Tj) and (122]) are positive for arbitrary k. From this, we conclude 
that the charged condensate background is stable w.r.t. small perturbations. All 
the group velocities obtained from f[2TT) and f[2~2~|) are subluminal. 

4 Brief discussions and outlook 

The spectrum of bosons exhibits a mass gap. As a result, contributions of the bosons 
into the specific heat of the substance at low temperatures would be suppressed as 
exp(— m 7 /T), or stronger. Since for the densities at hand m 7 ~ (30 — 90) KeV, 
(corresponding to ~ (3 — 9) • 10 8 K) the exponential suppression will be strong at 
lower temperatures. 

As the star cools further, the specific heat will be dominated by the contributions 
of the relativistic electrons near the Fermi surface, and, hence, would vanish with the 
temperature linearly, ~ T. This resembles the low temperature scaling in metals, 
but differs from that in insulating solids, or in the BE condensate of the neutral 
helium-4 atoms, where the low temperature specific heat is dominated by massless 
phonons, and scales as ~ T 3 . 

If charged condensate is realized in the helium dwarf stars, as we have argued in 
this work, there may be important consequences for the equation of state, cooling, 



2M 4 



mi 



+ k 2 



m A Q - 



rut 



2m 2 



(22) 
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rotation, as well as the asteroseismology of the helium dwarfs. These and related 
issues will be studied elsewhere. 
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